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Abstract 

It is shown that the effective inertial mass density of a dissipative 
fluid just after leaving the equilibrium, on a time scale of the order 
of relaxation time, reduces by a factor which depends on dissipative 
variables. Prospective applications of this result to cosmological and 
astrophysical scenarios are discussed. 


*Postal address: Apartado 80793, Caracas 1080A, Venezuela 
le-mail address: laherrera@telcel.net.ve 


1 



1 Introduction 


In classical dynamics the inertial mass is defined as the factor of proportion¬ 
ality between the three-force applied to a particle (a fluid element) and the 
resulting three-acceleration, according to Newton’s second law. 

In relativistic dynamics a similar relation only holds (in general) in the 
instantaneous rest frame (i.r.f.), since the three-acceleration and the force 
that causes it are not (in general) paralell, except in the i.r.f. (see for example 

i)- 

In this work we shall derive an expression for the effective inertial mass 
density of a dissipative relativistic fluid, which is valid just after the fluid 
leaves the equilibrium, on a time scale of the order of relaxation time. 

By “effective inertial mass” (e.i.m.) density we mean the factor of pro¬ 
portionality between the applied three-force density and the resulting proper 
acceleration (i.e., the three-acceleration measured in the i.r.f.). 

As we shall see, the obtained expression for the e.i.m. density contains a 
contribution from dissipative variables which reduces its value with respect to 
the non-dissipative situation. Such decreasing of e.i.m. density has already 
been brought out in the spherically symmetric self-gravitating case [^, the 
axially symmetric self-gravitating case and also for slowly rotating self- 
gravitating systems |^. Here we want to provide a general derivation for the 
e.i.m. density of a dissipative fluid, which is independent on the symmetry 
of the problem.On the other hand, in order to derive the expression for the 
e.i.m. density, it will be necessary to suppress gravitational contributions to 
four acceleration. This will be achieved by evaluating the equation of motion 
in a locally Minkowski frame, where local gravitational effects vanish. Then 
evaluating the resulting equation of motion in the i.r.f. we are led to an 
expression which takes the desired, ” Newtonian”, form. 

Force=e.i.m. x acceleration(proper). 

It is perhaps worth noticing that the concept of effective inertial mass is 
familiar in other branches of physics, thus for example the e.i.m. of an 
electron moving under a given force through a crystal, differs from the value 
corresponding to an electron moving under the same force in free space, and 
may even become negative (see p, ^). 

In order to illustrate our point and to establish the notation, we shall 
hrst consider, very briefly, the perfect fluid case. However before closing this 
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Section, the following comments are in order: 

It is worth noticing, that there is not a unique way to dehne the four 
velocity of a dissipative fluid. Indeed, we may take, among other possibilities, 
the four velocity vector as being related to the velocity of energy transport 
(Landau-Lifshitz ||^) or to the velocity of particles (fluid element) as done by 
Eckart (0), both choices are of course perfectly equivalent, and both lead to 
different dehnitions of rest frame. In the former case there is not energy flux 
in the rest frame (the components of the energy tensor vanish), whereas 
in the later the spatial components of the particle current four vector vanish 
in the comoving frame. In this work we adopt the Eckart choice. 

This is motivated by the fact that the very concept of e.i.m. as dehned 
above, requires (and makes sense only for) such choice (Eckart). Indeed, 
observe that the concept of e.i.m. requires the existence of a frame where 
the three-force applied to a particle and the resulting three acceleration are 
parallel, and this only happens in the i.r.f. where a given fluid element is at 
rest. But this is the Eckart frame by dehnition. 

In the Landau-Lifshitz frame we cannot say that such relation between the 
applied three-force and the resulting three accelerations follows, and therefore 
in such case the very meaning of e.i.m. is uncertain. 

2 The e.i.m. density of a relativistic perfect 
fluid 

Let us consider a relativistic perfect fluid, whose energy-momentum tensor 
takes the usual form 

= {p + p)u^u'^ - pg^'^ (1) 

(we use relativistic units and signature —2) 

Then the equation of motion yields 

(p + p)a“ = (2) 

with 

^ ( 3 ) 

u>^ = a>^ = ( 4 ) 
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where colon and semicolon denote partial and covariant derivatives, respec¬ 
tively. 

Observe that for the sake of generality we are including gravitational held, 
although, as already noticed, the hnal expression for e.i.m. density will be 
independent of gravity. 

Eq.(^ suggests that {p+p) plays the role of the e.i.m. density of the huid 
(which is in fact the case). However, we recall that the applied three-force 
and the resulting three-acceleration are paralell only in the i.r.f., therefore, 
in order to conclude that (p -|- p) represents the e.i.m. density, (@) has to 
be evaluated in such a frame. Furthemore, it should be observed that a" 
contains both, gravitational (in the self-gravitating case) and “kinematical” 
contributions. Since in the concept of e.i.m. only the later should appear, we 
have to evaluate (^) not only in the i.r.f., but also in the locally Minkowskian 
frame (l.m.f.) where local gravitational effects vanish. Formally, this may be 
achieved, either by introducing locally Minkowskian coordinates or equiva¬ 
lently, by considering a tetrad held attached to such l.m.f.. Denoting by a 
tilde the components of tensors in the l.m.f.,equation (0) reads 

{p + p)a^ = h^"''p^u (5) 

Evaluating (^) in the i.r.f., it becomes clear that {p + p) dehnes the e.i.m. 
density of any huid element (a well known fact, see for example |^). 

In order to better illustrate this point, let us consider the spherically 
symmetric case. 

Thus, let the line element be 

ds^ = e'"dt^ — e^dr"^ — r‘^{d9‘^ + sin^ dd(f)^) (6) 


with V = iy{r,t) and A = A(r, t), then locally Minkowski coordinates may be 
introduced by |T^ 


dT = e^^'^dt dx = e^^'^dr dy = rd9 dz = rsin9d(j) 
satisfying 

ds^ = dT'^ — dx^ — dy‘^ — dz^ 

Next, the four-velocity vector for any huid element, in the frame of (^ 


= 


\/l — ’ \/l — 


, 0,0 


(7) 

( 8 ) 

(9) 
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where to is the radial velocity of the fluid element as measured in the l.m.f. 
Therefore, in the l.m.f. the components of are 


= 


5t^ 




+ 


uj5^ 




\/l — 


LO^ 


( 10 ) 


and the components of the four-acceleration in the l.m.f. read 


Then, feeding back ([TT|) into (H), and evaluating in the i.r.f. {u = 0), we 
obtain the “Newtonian-type” equation 


{p + p)uj^T =-P,x (12) 

indicating that {p + p) represents the e.i.m. density of the fluid. 

Let us now turn to the dissipative case. 


3 The e.i.m. density of a dissipative fluid 

For simplicity we shall consider the case of pure heat conduction, neglecting 
viscosity terms. Then, the energy-momentum tensor takes the usual form 

= (p + p)u^u^ - + QtiUu + (13) 

and the equation of motion, projected onto the plane orthogonal to reads 


(p + p)a“ - = 0 


(14) 


where g^, satisfying = 0, denotes the heat flow vector. 

We have now to adopt a transport equation for g^. Although this is still 
a matter of discussion, it should be borne in mind that some kind of hyper¬ 
bolic dissipative theory should be applied in order to avoid violation of 
causality and other undesirable consequences derived from the (parabolic) 
Eckart-Landau approach. Furthermore, since we are interested in the tran¬ 
sient regime between two equilibrium states we are compelled to use a hy¬ 
perbolic theory, for it is known that the parabolic theory assumes that the 
system is relaxed at all times. 
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Thus, we shall consider the Israel-Stewart transport equation 
though any hyperbolic equation yielding a Cattaneo-type eqi 
the non-relativistic limit would lead to the same result) 

1 f 'T~i y ^ \ 

- Ta,) - -kT^ j 


IT| (al¬ 


ii equation in 


(15) 


where r, k, T and denote the relaxation time, the thermal conductivity, 
the temperature and the vorticity tensor, respectively. 

Then, feeding back (^) into (|T^ one obtains 

kT 


/ K \ 1 

[P + P-—)- P,-h^" 


2r 


-kT^ 


TW 


+ g“0 + gV“ = 0 


;/3 


(16) 


with 


Q = <, 


as usual denotes the shear tensor, and where the relation 




(17) 


has been used. 

We shall now consider that our system is initially, either in hydrostatic 
equilibrium without dissipation, or very close to equilibrium with “small” 
dissipation (quasi-st at ionary regime). Next, let us suppose that at some 
moment (say t = t) it leaves the initial state entering into a dissipative 
regime. 

Then, “immediately” after departure from equilibrium (or quasi-equilibrium), 
where “immediately” means on a time scale of the order of r (or smaller), 
the last four terms in (|I§ can be neglected, and we have 

/ k' !' \ 1 

a^{p + p- — j - = 0 (18) 

The components of ([T^) in the l.m.f. are 

(p + p - = 0 (19) 
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However, in the l.m.f., the components of the heat flux vector in the i.r.f., 
in the quasi-stationary regime, satisfy the “Maxwell-Fourier” equation 




( 20 ) 


Indeed, just after leaving the equilibrium, the absolute value of the heat 
flow vector is the same as that corresponding to the quasistationary regime, 
and furthermore in that regime (in the l.m.f. and in the i.r.f.) both the 
kinematical and gravitational contributions to the four acceleration van¬ 
ish, thereby justifying (p0|). Therefore , just after leaving the quasi-stationary 
regime, the last term in (|^ vanishes, and this equation becomes 


[p + V- 

implying that the e.i.m. density is given by 

kT 

P + P - 

r 


or, 


with 


{p + p){l-a) 

kT 


a = 


t(p + p) 


( 21 ) 

( 22 ) 

(23) 

(24) 


4 Discussion 

We have seen that just after leaving the equilibrium the e.i.m. density of 
the dissipative fluid reduces its value (as compared to the non-dissipative 
situation) by the factor 

1 - a (25) 

giving rise, in principle, to the possibility of vanishing {a = 1) or even neg¬ 
ative (q; > 1) e.i.m. density. As mentioned in the Introduction, such effect 
was already found in some specihc examples, here it appears as a general 
result for any dissipative fluid (self-gravitating or not). 

It should be noticed that since the system is evaluated immediately (in 
the sense explained above) after leaving the equlibrium (or quasi-equilibrium) 
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then the physical meaning of thermodynamical variables (their values), are 
extrapolated from the state of equilibrium, and therefore any discussion 
about difficulties of interpreting such variables in a transient regime, is out 
of the point. 

Also, it should be observed that causality and stability conditions hin¬ 
dering the system to attain condition a = 1, are obtained on the basis of a 
linear approximation, whose validity, close to the critical point {a = 1), is 
questionable 
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At any rate, examples of fluids attaining the critical point and exhibiting 
reasonble physical properties have been presented elsewhere [|^ pA[| . 

In order to evaluate a, let us turn back to c.g.s. units. Then, assuming 
for simplicity p + p k. 2p, we obtain 


kT 


_ m 

t{p + p) ^ [t][p] 


X 10 


-42 


(26) 


where [«:], [T], [r], [p] denote the numerical values of these quantities in 
erg. cm~^ K, s and g. cm~^, respectively. 

Obviously, this will be a very small quantity (compared to 1), unless 
conditions for extremely high values of k and T are attained. In this respect. 


observe that although small values of r increases (^61) , they are of little 
help since ( p^ is valid only on a time scale of the order of r. Therefore, a 
decreasing of e.i.m. density with physically relevant consequences requires 
values of (@) close to the unit, due to large values of k and T but non- 
negligible values of r. 

At present we may speculate that a may decrease substantially (for non- 
negligible values of r) in a pre-supernovae event 

Indeed, at the last stages of massive star evolution, the decreasing of 
the opacity of the fluid, from very high values preventing the propagation 
of photons and neutrinos (trapping [|T^]), to smaller values, gives rise to 
radiative heat conduction. Under these conditions both n and T could be 
sufficiently large as to imply a substantial increase of a. Indeed, the values 
suggested in |]^ {[k] ^ 10^’^; [T] 10^^; [r] 10“^; [p] ^ 10^^ ) lead to 

ct 1. The obvious consequence of which would be to enhance the efficiency 
of whatever expansion mechanism, of the central core, at place, because of 
the decreasing of its e.i.m. density. At this point it is worth noticing that the 
relevance of relaxational effects on gravitational collapse has been recently 
exhibited and stressed (see 0, and references therein) 














It is also worth noticing that the inflationary equation of state (in the 
perfect fluid case) p + p = 0, is, as far as the equation of motion is concerned, 
equivalent to a = 1 in the dissipative case (both imply the vanishing of the 
e.i.m. density). 

Also observe, that if we impose the equation of state p + p = 0, in the 
dissipative case, then the resulting e.i.m. density is negative (-^), implying 
that the attractive gravitational force leads to an overall expansion. 

Finally, it is worth stressing that it is the first term on the left in ([T5|) , 
the direct responsible for the decreasing in the e.i.m density, all terms includ¬ 
ing kinematical variables vanishing within the time scale considered (after 
leaving equilibrium). Therefore any hyperbolic dissipative theory yielding 
a Cattaneo-type equation in the non-relativistic limit, is expected to give a 
result similar to the one obtained here. 
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